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I. INTRODUCTION 

Decoherence through environmental coupling is a 
crucial ingredient for understanding the quantum-to- 
classical transition Em. As a prototypical problem, one 
may consider the motion of a quantum particle interact¬ 
ing with an ideal gas environment. The interaction with 
the environment selects a preferred set of states, the so- 
called pointer states M- These special states of the 
particle remain pure for a relatively long time, whereas 
their superpositions decohere very fast into a mixture. 
For simple master equations the pointer states are known 
to be localized wave packets moving in phase space. 
In the semiclassical limit they turn into arbitrarily sharp 
peaks moving along classical phase space trajectories. As 
such, pointer states are ideally suited to explain the emer¬ 
gence of classicality within a quantum framework. 

The name pointer states derives from the fact that 
their characteristics matches the behavior of a pointer 
of a measuring apparatus. Namely, if a quantum system 
is in an eigenstate of the observable we measure with an 
apparatus its pointer will turn to a specific value on the 
dial and stay there. On the other hand, if the system is in 
a superposition of two nondegenerate eigenstates one will 
not find the pointer of the apparatus to be in a superpo¬ 
sition of the two corresponding values, as prescribed by 
unitarity, but pointing to either one or the other position 
with probabilities determined by Born’s rule. 

In this article, we identify the pointer states of the mas¬ 
ter equation of quantum Brownian motion (QBM) and 
derive their equations of motion. Specifically, we under¬ 
stand as QBM the dissipative dynamics of a free particle 
resulting from many small environmental collisions as de¬ 
scribed by a Lindblad type master equation E3:EI|. One 
possibility to obtain this master equation is to consider 
a particle linearly coupled to a bath of harmonic oscil¬ 
lators E2“E1| and to take the high-temperature limit of 
the reduced particle dynamics; one arrives at the master 
equation of the Caldeira-Leggett model (1,111 Eg. To 
get a Lindblad master equation an additional term has 
to be appended EMI making the equation completely 
positive. Alternatively, the Lindblad master equation of 


QBM can be obtained by considering a dissipative limit 
of the general particle motion in a gaseous environment 
as governed by the so-called quantum linear Boltzmann 
equation (TK^ . 

By considering quantum trajectories using a particular 
stochastic process in Hilbert space, the so-called pointer 
state unraveling, we find the pointer states of QBM to 
be rotated Gaussians in phase space, in line with the 
encounter of Gaussian pointer states in other models of 
linear environmental coupling d d . Moreover, we show 
how superpositions of separated wave packets decohere 
into mixtures of pointer states with weights given by the 
absolute value of the initial amplitude in the superposi¬ 
tion (i.e. Born’s rule). Finally, we derive the stochastic 
equations of motion for the location of the pointer states 
in phase space and we show that they turn into the equa¬ 
tion of classical diffusion in the semiclassical limit, thus 
providing a comprehensive picture of the quantum-to- 
classical transition for the important dissipative process 
of QBM. For times of evolution much greater than the 
characteristic decoherence time, it is then justified to ap¬ 
proximate the dynamics by a mixture of pointer states 
moving on trajectories in phase space, instead of calcu¬ 
lating the whole evolution of the density matrix. 

Our method of using the pointer-state unraveling to 
identify the pointer basis and to provide a perspective 
on the quantum-to-classical transition for particle motion 
has already been demonstrated for the non-dissipative 
case of pure collisional decoherence d, m, HI. In that 
case, the pointer states are exponentially localized, rather 
than Gaussian, and the emerging classical equations of 
motion are Newton’s equations. In the present case of 
dissipative dynamics we derive a stochastic equation of 
motion for the pointer states in order to get the Langevin 
equation corresponding to classical diffusion. This re¬ 
quires a subtle analysis of the dynamics of the quantum 
trajectories. 

The paper is organized as follows: At first, in sec¬ 
tion im we give a general definition of pointer states. In 
the following section Hill we introduce and briefly discuss 
the Lindblad master equation of QBM and introduce its 
pointer state unraveling. The careful examination of this 
particular quantum trajectory approach is the topic of 
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the remainder of the paper: The deterministic part of 
the unraveling is analyzed in section HVl providing us with 
the Gaussian pointer states of QBM as well as with the 
deterministic part of the classical equations of motion. 
By looking at the full unraveling, i.e. the deterministic 
and stochastic part, we show in section^that any super¬ 
position of separated wave packets turns into a mixture 
of pointer states with the expected probabilities. In ad¬ 
dition, we identify the trajectories of the pointer states; 
they exhibit a stochastic motion in phase space and yield 
diffusion in the semiclassical limit. Section EH contains 
our conclusions. 


II. POINTER STATES 


equation associated with the master equation. This non¬ 
linear pointer state equation (NLPSE) can be obtained 
heuristically by identifying the pure state evolution clos¬ 
est to the master equation evolution (^. [^. 

The pointer state unraveling, a piecewise deterministic 
stochastic process in Hilbert space, connects the solution 
of the master equation with the NLPSE, since the latter 
describes the deterministic part of the unraveling. In the 
next section, we present the master equation of QBM 
together with its pointer state unraveling. 

III. QUANTUM BROWNIAN MOTION 
A. Master equation of quantum Brownian motion 


The evolution of a Markovian open quantum sys¬ 
tem can be described by a Lindblad master equation 
[laiill,!!!. For this equation to exhibit pointer states 
we expect it to display a time scale separation into a 
short decoherence time tdec and a much longer dynami¬ 
cal time scale which has also a classical meaning. The set 
of pointer states is characterized by the property that su¬ 
perpositions of mutually orthogonal pointer states decay 
into a mixture on the decoherence time scale whereas a 
single pointer state remains relatively stable and changes 
on the longer time scale only. 

Specifically, for a superposition of pointer states 

N 

IV'o) = (1) 

a=l 

the solution of the master equation for times much 
greater than tdec is well approximated by the mixture 
of pointer states 


The master equation of QBM describes the quantum 
state of a massive marker particle linearly coupled to a 
high-temperature gas environment and experiencing de¬ 
coherence as well as dissipation. For the sake of sim¬ 
plicity, we restrict ourselves here to the one dimensional 
case and use dimensionless variables by introducing time, 
length and momentum scales 

T=-;^, ^ = P = ^mkBTenv (4) 

They involve the mass of the particle to, the tempera¬ 
ture of the surrounding environment Tenv, and the fric¬ 
tion constant 7 . These scales do not depend on h and are 
therefore convinient in describing the semiclassical limit 
of QBM. In these units the canonical commutator takes 
the form [x, p] = ij k and the momentum operator in po¬ 
sition representation is (a;|p|<I>) = —i/K9x4>(x), where we 
introduced the dimensionless parameter 


N 

Po = |V'o)(V'o| Pt-'^ Prob(Q!|po)ka(t))(7ra(t)|. 

a—1 

( 2 ) 

A proper set of pointer states forms a basis, is indepen¬ 
dent of the initial state, and depends only on the envi¬ 
ronment and the coupling operators. As a consequence 
and according to Born’s rule, the only dependence of the 
mixture m on the initial state po is via the probabilities 

Prob(a|po) = Tr (po| 7 ra( 0 ))( 7 ra( 0 )|) = |ca(0)k (3) 

The time dependence of the pointer states \ijja{t)) occurs 
on the longer time scale only. Once the pointer states and 
their time dependence are characterized, the dynamics of 
the master equation is fully captured by Eq. ([2]) for times 
greater than tdec- 

There are different possibilities to identify the pointer 
basis. A well-known one requires that the pointer states 
should be the least entropy producing states of the sys¬ 
tem; this is the so called predictability sieve [ 2 ^, [?7| . 
Here, we follow another approach, in which the pointer 
states are identified as soliton-like solutions of a nonlinear 


K = 


kBTe, 


2'yh 


(5) 


The semiclassical limit is then obtained by sending k —>■ 
00 . 

The Lindblad master equation of QBM is obtained by 
using a single Lindblad operator 

L = \/2 -I- ( 6 ) 

and the Hamiltonian 

H = (y + ^ {x, p}^ . (7) 


The anti-commutator {x, p} /4 describes a rescaling of 
the particle energy due to the environmental coupling. 
Inserting and ([7]) into the general Lindblad master 
equation flO. 



1 

\h . 


H,p 


LpL'l' - 



( 8 ) 
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we arrive at the Lindblad master equation of QBM (T^ . 

m 




-if [P^^] - if [><:{P>P}] - 
- ^[P>[P:P]]- 


(9) 


Equation © is the Lindblad generalization of the fa¬ 
mous Caldeira-Leggett master equation [llllllll, hav¬ 
ing in common the first three terms on the right-hand 
side. They describe free motion, friction, and monien- 
tnm diffusion, respectively. The last term of ([5]), which 
ensures complete positivity, can be identified with a posi¬ 
tion diffusion. This position diffusion term is not present 
in the Caldeira-Leggett master equation, and its influ¬ 
ence decreases with growing n. 

When writing © in the Wigner phase space represen¬ 
tation one can draw the semiclassical limit, k —>■ oo, to 
see the formal analogy between the evolution equation 
for the Wigner function and the Fokker-Planck equation 
for the classical phase space distribution of classical dif¬ 
fusion Bill. Further analogies to the classical case of 
Brownian motion can be seen if one considers the mean 
values and variances of the position and momentum op¬ 
erators [ 2 ^. Most prominently, one finds that the mean 
of the momentum operator is exponentially damped and 
that the variance of the position operator grows linearly 
for asymptotically large times; both properties are char¬ 
acteristics of classical Brownian motion. 


B. Unraveling of the master equation 


In order to derive the pointer states of QBM and 
their motion in phase space, we replace the deterministic 
master equation evolution for the density matrix ([HI) by 
an equivalent stochastic pure state evolution in Hilbert 
space, the pointer state unraveling. The equivalence of 
an unraveling to a master equation is provided by the 
fact that the ensemble average over the quantum traje- 
tories yields the solution of the master equation [^, [ 28 l - 
[33|. Since the decomposition of a density matrix into 
a mixture of pure states is not unique, many ensembles 
of quantum trajectories will recover the same solution of 
the master equation, hence the choice for an unraveling 
is ambiguous. One distinguishes between continous un- 
ravelings based on a Wiener process , where 

quantum state diffusion is the most prominent, and piece- 
wise deterministic unravelings [1, |J[2^, [2^, [sj, [s^ based 
on a Poisson process. 

We will use here a piecewise deterministic unraveling 
which was already employed in the case of collisional de¬ 
coherence to determine the pointer states and their dy¬ 
namics ^ [^, [3^ . In our case of QBM the deterministic 


part is formed by the nonlinear operator 
(p^-p^) 

+ (k - (x - xf'j + ^ “ (p “ 

+ ^ (p + p) ^ IV'), 

( 10 ) 

with expectation values denoted as H = (A) = (^/;|A|i/)), 

variances Va = A^—A and the covariance Cxp = {x,p}— 
2xp. Due to the time dependence of the state \ijj), all 
expectation values are time dependent. The time variable 
is omitted here for the sake of readability. 

The evolution equation 

^IV') = N[V']|^/>), (11) 

is called the nonlinear pointer state equation (NLPSE) 
because it is expected to exhibit soliton-like solutions 
which we will identify as the pointer states. Soliton-like is 
understood here in the sense that the shape of the wave 
function is time-independent in phase space. As men¬ 
tioned above, the NLPSE m can also be obtained by 
looking for the particular pure state equation that mimics 
the master equation most closely [y, . 

To obtain the whole unraveling a stochastic jump part 
must be added to the deterministic NLPSE m yielding 

Id^;) = NM|^)df-h [ - 1 ] l^)dA^, (12) 

with the nonlinear jump operator 

J[V'] = \/2 I^k(x-x)- b ^ (p-p)^ . (13) 

The stochastic Poisson increment dA^ can take values 
dN G {0,1} and satisfies the Poisson increment rule 


dN^ = dN. (14) 

The nonlinearity of the jump operator (1131) arrises be¬ 
cause of its dependence on the position and momentum 
expectation values which depend on the state |'(/')• To 
fully define the above Poisson process we have to specify 
the ensemble average of the Poisson increment. It is de¬ 
termined by the rate r = (^/>|jl [■!/)] J[i/)]|'0) of the process, 
given in our case by 


E[dN] 

dt 


= r = 2 k^Vx + — 


2 ’ 


with E [•] denoting the ensemble average. 


(15) 
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A straightforward calculation shows that this unrav¬ 
eling indeed recovers the master equation (jH]) in the en¬ 
semble average It is also straightforward to verify 
that the jump rate in this unraveling is proportional to 
the change of purity of the state 

(16) 

This finding fits with the expected behavior of a pointer 
state unraveling: Once the pointer state has been reached 
due to the deterministic part of the unraveling (as deter¬ 
mined by the NLPSE) it should be relatively stable. This 
implies that the change of purity of the pointer state is 
small and so is therefore the total jump rate (fTOl) . As a 
consequence, once a pointer state is formed in a quantum 
trajectory of this unraveling, it remains a pointer state 
as long as possible. 


IV. POINTER STATES OF QUANTUM 
BROWNIAN MOTION 

In this section we present the pointer states of QBM 
as soliton-like solutions of the NLPSE (ED. After that 
we analyze the action of the NLPSE on a superposition 
of separated wave packets, which will help us to obtain 
Born’s rule in section El 


A. Single wave packets 

A lengthy but straightforward calculation confirms 
that the NLPSE dill) is solved by a Gaussian state of 
the form 


ip{x) 


1 

(2^14)'/^ 



{x — x)'^ 

414 


(1 


iMiCxp) 


+ in {x 


x)p + i4> 


(17) 


where is a time dependent global phase. This requires 
that the first moments and the variances satisfy the fol¬ 
lowing closed set of differential equations. 


d_ 

d_ 

-14 = Cxp - K (4k:"14 - 1) + ^ ~ 


dt 


16k^ 


_d_ „ 1 + 8 — Cx 

dt^ 


^Cxr, = 






(18) 

(19) 

( 20 ) 
( 21 ) 


To see this, one inserts the ansatz ([T7| into the NLPSE 
ED and compares coefficients in powers of a;. A differ¬ 
ential equation for I 4 is readily derived by noting that 


the variances of the Gaussian state ED are related via 
4VxVp = ^ + ( 22 ) 

Kj 

The equations of motion (fT^ and (|19ll for the first mo¬ 
ments exhibit an exponential momentum damping, which 
reflects the friction experienced by the Brownian particle 
through the interaction with the gas environment. This 
movement does not affect the form of the wave packet 
but only its position in phase space. The equations of 
motion (EHD and ED for the variances, however, drive 
the width and orientation of the wave packet in phase 
space towards the stable fixed values 


14 ,ps - 


C. 


aip,ps 


1 -f + (64k2 -H 1) V16k2 +J 
8kV16k2 -h 1 
4k — \/ -1-1 — 1 
k\/16k^ -I- 1 


(23) 

(24) 


The label “ps” denotes the pointer state, and 14,ps can be 
calculated with the help of (1^^ . These fixed values are 
obtained by solving the nonlinear equations which arise 
by setting dl4/dt = dCxp/dt = 0; their stability can be 
confirmed by a linear stability analysis of (1201) and ED 
around 14 ,ps and Cxp,ps- 

A wave packet with the fixed values 14,ps and l^.ps 
is a soliton-like solution in the sense that the form of 
the wave function in phase space does not change with 
time. We will call this solution a pointer state of quantum 
Brownian motion and denote it by [tt {x,p)), where x and 
p parametrize the position and momentum coordinates 
of the pointer state. As Gaussian wave functions they 
form an overcomplete basis set in Hilbert space. From 
Eqs. (I22l) - ((24l) one deduces the asymptotic widths of the 
pointer states as k —>■ 00 , 


14 


tc,ps 


14,ps 


Cx 


1 

24^’ 

1 

1 

K 


(25) 


The pointer states thus get more and more localized in 
phase space as one goes deeper into the semiclassical 
regime (fi —>■ 0) characterized by k —)> 00 , see Eq. ([S]). 
This is illustrated in Fig. [TJ 


B. Superposition of wave packets 

We now analyze the action of the NLPSE ED on a 
superposition of separated wave packets. We will see that 
the NLPSE suppresses such superpositions until only a 
single localized wave packet remains, which eventually 
turns into a pointer state \tt {x,p)) introduced before. 
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Figure 1. Pointer state variances 14,ps in position (solid 
line) and 14 ,ps momentum space (dashed line) and covari¬ 
ance Cxp,pB (dotted line) as a function of k. For k 3> 1, i.e. in 
the semiclassical limit, all curves follow a potential law. The 
scales are defined in Eq. 


We consider an initial superposition state of the form 
\i’) = (26) 

3 

where the wave packets |'0j) are well separated in phase 
space, orthogonal, and normalized with weights 

= |c,(0)|% (27) 

summing up to unity. We drop the time argument of the 
weights and amplitudes for better readability, and we de¬ 
note expectation values of the wave packet components 
as Aj = |A|'0j), and the variances and covariance ac¬ 

cordingly. 

The are required to be well separated wave packets 
in the sense that the orthogonality condition 

(V'i IAI IA1) Vj, fc, (28) 

holds for the operators A g{I, x,x^, p, p^,xp -|- px}. The 
shape of the \il3j) can be arbitrary, but their widths in 
phase space shall be close to that of a pointer state. In 
particular, we require that the variances 14,j and Vpj of 
the wave packet components \il}j) exhibit the same depen¬ 
dence (1^51) on K as the pointer state. The orthogonality 
condition (E51) can then be fulfilled arbitrarily well by 
increasing the parameter k, which decreases the pointer 
state width in phase space according to (1^51) . and thus 
further orthogonalizes the 

For the following, it is useful to express the expectation 
values of the superposition state I'i/’) with the help of 


those of the constituting wave packets |'0j): 

^=J2wjXj, 

j 

p = j2w,p^, 
j 

3 3,k 

Vp=J2 ^3^P,3 (Pj -Pkf ^ 

3 3,k 

C^P = {Xj - Xk) {Pj -Pk) ■ 

3 hfe 

(29) 

By inserting the superposition state (I^Sl) into the NLPSE 
m we get the expression 

3 3 

(30) 

where the nonlinear operator on the right-hand side con¬ 
tains expectation values with respect to |'(/')■ From this 
equation we deduce a set of coupled equations for the 
normalized states |'0j ) and their weights Wj = |cjp: 

(x - Xj) + ^ (p - Pj)^ 

X (^k{x-Xj) - ^{p- pj)^ (31) 



+ g i^p ^p,3 (p Pj) ) ^ • (32) 

Note, that the nonlinear operator appearing in pil) now 
only contains expectation values with respect to I'll’j)- 
Moreover, if there is only one constituent in the “super¬ 
position”, Eqs. (I5T]) and (1321) turn consistently into the 
NLPSE (ITT]) and the trivial evolution dwj /dt = 0 for the 
weight. This is seen most clearly by looking at the dif¬ 
ferences of the expectation values, which vanish in that 
case. 

The coupled weight equations (15^ exhibit a fixed point 
if a single wj is equal to one, while all others are zero. 
It can be checked by a linear stability analysis that this 
fixed point is a stable one. Thus, Eq. (|5^ expresses the 
fact that the NLPSE suppresses superpositions of pointer 
states, since the dynamics of the weights always ends up 
in the stable fixed point. 

To illustrate the suppression of a superposition state, 
we consider a superposition of two pointer states. The 
weight of the first component I'i/’i) then follows the dif- 
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ferential equation 


dt 


wi = {2wi — 1 ) W 1 W 2 ( 2k'^ {x 2 — xi)^ + - {p: 



We see immediately that wi gets greater if and only if it is 
already the greater of the two weights, leading eventually 
to the suppression of the other component. If there exists 
a superposition of many wave packets, which, in addition, 
all have different widths the situation becomes far more 
complex and can not be captured intuitively. However, 
the eventual decay of the superposition into one of the 
components is certain, though it is not easily predictable 
which component of the superposition will survive in the 
course of the evolution due to Eq. (15^ . 


V. UNRAVELING QUANTUM BROWNIAN 
MOTION 

In this section we return to the unraveling of QBM, 
which consists of the NLPSE as deterministic part, com¬ 
plemented by stochastic jumps. This stochastic part is 
neccessary to produce the required ensemble of quantum 
trajetories since the NLPSE alone would always lead to 
the same asymptotic pointer state for one particular ini¬ 
tial state. We first show how the jump process leads to 
stochastically selected asymptotic pointer states, which 
occur with relative frequencies according to Born’s rule. 
After that, in section IVB[ we discuss the action of the 
unraveling on a single pointer state. The NLPSE de¬ 
scribes pointer state trajectories which are exponentially 
damped in momentum. We will see how the stochastic 
part induces a random walk of the pointer state trajec¬ 
tory displaying the expected diffusive behavior. In the 
semiclassical limit this will eventually yield classical dif¬ 
fusion. 


A. Stochastic dynamics of a superposition 

Let us consider the action of a single jump on the 
superposition state dMl) of section lIVBi consisting of 
wave packets lijjj) that are well separated in phase space, 
according to the orthogonality relation (1^51) . and with 
widths comparable to the pointer state width, 

w')=^=i: <=;«>■ (34) 

Here, we introduced the normalization factor (J^J) = 
J[V’]|V’) s-nd the normalized wave packets 






( 35 ) 


involving the normalization factor 

(Vx,j + {x-Xjf^ 

(36) 

I 12 

In the following, only the squared moduli re' = |c'| of 
the new coefficients will be required. One readily finds 



(37) 


Note that the new wave packet components \tjjj) can 
be safely assumed to be still separated and localized be¬ 
cause the jump operator modifies the shape only linearly 
in X and p. The superposition state (1341) after the jump is 
therefore again a superposition of separated wave pack¬ 
ets, but with different weights. The effect of a jump can 
thus be approximately accounted for by a reshuffling of 
a finite number of weights Wj in the superposition. We 
confirmed this by a numerical implementation of the un¬ 
raveling (fl^ by combining the Crank-Nicolson method 
with the split-operator technique (dSl - ld^ . 

From the action of the NLPSE on a superposition state 
we calculated in section HVBI the deterministic evolution 
(1321) of the weights Wj. Knowing the new weights dsn 
due to a jump, we are now in a position to construct 
a stochastic process for the weights of the wave pack¬ 
ets where the jumps occur with the same rate as in the 
original unraveling: 


dwj = Wj ((j^J) - (jl'J)j) dt 

-f W4 ( - 1 I dW (38) 

Here, we have written the deterministic part (I32|) more 
conveniently with the help of the jump rate r = (J u) (fT51) 
and the normalization factor (J^J)j (|36l) . The Poisson 
increment diV has the same ensemble average E[dV] = 
r dt as the pointer state unraveling m- 

The set of equations (p8)) represent the evolution of the 
weights in a quantum trajectory. Each quantum trajec¬ 
tory ends up asymptotically in a pointer state, when all 
weights but a single one vanish. From this we deduce 
that the relative frequency of finding a particular pointer 
state \tt (xjjPj')) is equal to the probability of the as¬ 
sociated asymptotic state in the ensemble of quantum 
trajectories. This latter probability is evaluated easily 
by carrying out the ensemble average of (1551) . 


E [dwj] = 0. 


(39) 


This implies E [wj(t)] = const, and, in particular. 


E [wj{t ^ 00 )] = E [wj{0)] = |cj(0)p, 


(40) 
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where Cj (0) are the initial amplitudes of the wave packets 
in the superposition state Equation (HOI) confirms 

Born’s rule ([H), as expected for a superposition of sepa¬ 
rated wave packets. 


B. Stochastic dynamics of a single wave packet 


So far, we investigated the dynamics of the unravel¬ 
ing if the initial state is a superposition of localized wave 
packets. We confirmed in the previous section that the 
probability to end up asymptotically in a certain pointer 
state is given by its weight in the initial superposition 
(Born’s rule). This process takes place on the fast de¬ 
coherence time scale tdec after which the system is well 
described by a pointer state. In terms of the pointer state 
unraveling, this point is reached when each quantum tra¬ 
jectory in the ensemble has turned into its asymptotic 
pointer state. 

Since the underlying quantum dynamics is diffusive 
[ 2 EI it is natural to describe the statistics of the posi¬ 
tion X and the momentum p of the pointer state by the 
phase space diffusion 




Bi2 W dWi \ 
B 22 ) \dW 2 


(41) 


with the Wiener increments dWi and dW 2 obeying the 
rules 


deterministic motion is interrupted by a jump, described 
by the jump operator (1131) acting on the pointer state 

dm), 

-tfj'ix) = (x|J[7r(x,p)]|7r(s,p)) 



1 — IkCx 
8k^V, 


a;,ps 


(x 


x) 


(x|7r(x,p)). 


(44) 


Here, 


r — 

' ps — * ai,ps 


2 


(45) 


is the jump rate m of the pointer state and V^^ps 
and lp,ps are the widths of the pointer states defined 
in Eqs. (H^ - (IM1) . The resulting state is a symmetric 
double-peaked wave packet. This symmetry reflects the 
assumed symmetry of the state before the jump. Since in 
reality the system state is never exactly symmetric even 
small asymmetries will trigger the emergence of a single 
wave packet under the NLPSE, see Fig. [5] 

Which of the two subpeaks survives depends on the 
asymmetry before the jump. We assume equal probabili¬ 
ties and that the single wave packet is restored sufficiently 
fast by the NLPSE. This results in an effective jump of 
the pointer state’s first moments in phase space by the 
distances 


E[dW] = 0, 
dWidWj = 5ijdt. 


(42) 


Equation (|4ip is a special case of the two-dimensional 
Ornstein-Uhlenbeck process with constant coefficients 
(iH and the derivation of its ensemble variances and co- 
variance is presented in the appendix. The different dif¬ 
fusion constants Dx, Dp, and Dxp describing position, 
momentum, and covariance diffusion are related to the 
Bij by 


Dx — Bfi + Bi2, 

Dp = Bl, + BI 2 , (43) 

Dxp = B11B21 + Bi 2 B 22 - 

We proceed to show how the stochastic jumps acting 
on the motion of the pointer state lead to phase space tra¬ 
jectories which turn into classical diffusion in the semi- 
classical limit K ^ 00 . To achieve this, we derive a simple 
model allowing the description of the pointer state trajec¬ 
tories by the phase space diffusion (HD) and then compare 
it to a numerical study. 


jx — \x x| — \/2Vx,ps, 
jp = \p' -P\ = \/2Pp.ps. 

The widths jx and jp are calculated from the positions 
of the two peaks in (H^ . Since the direction of the jump 
is positive or negative with equal probability, we end up 
with the two possible jumps in phase space 


each of them characterized by a Poisson process dA^^ with 
jump rate 


_ E[djVfc] _ ^ 
dt ~ 2 ' 


(48) 


We are now in the position to write down a stochastic 
differential equation for x and p, 




2 

dt + '^jkdNk. 

/c=l 


(49) 


1. Analytic model 

Consider a pointer state moving on the phase space 
trajectory governed by the NLPSE (see Sec. IIV Al) . This 


It consists of the deterministic evolution (1181) and (IT^ de¬ 
rived in section HV Al and the stochastic jumps (H71) with 
corresponding jump rates (1481) . By using the Poisson rule 
dNkdNi = SkidNk and Eq. (H51) . one can easily derive all 













Figure 2. Simplified jump process: A jump turns tfie wave packet at position x (a) into a slightly asymmetric double-peaked 
structure (b). The NLPSE then suppresses one of the peaks resulting in a single wave packet at the position x' of one of the 
double-peaks (c). A similar picture applies in momentum space. 


moments of the stochastic process (H^ . 


E 

\dxf"-^' 

= 0, 

E 

\dxf'^' 

E 

(dp)""-'' 

= 0, 

E 

'(dp)"" 



E[(da:dp)”] = (jx jp)" J-psdt, 


with n > 1. If one is interested in calculating moments 
up to second order only, the stochastic jump process 
can be well approximated by the phase space diffusion 
(mu upon identifying the diffusion constants (1431) with 

^00 — jx'^psj 

Dp ~ ( 51 ) 

Dxp — jxjp'^ps- 

Inserting the jump widths (HSl) and the jump rate (H51) 
into Eqs. ED) one gets their dependence on k as k —>■ cxd 


Dp - 2, (52) 

1 

D-P - 

We see that both and D^p vanish in the semiclassical 
limit K —>■ oo so that only the momentum diffusion Dp = 
2 contributes and the stochastic process (HD) thus reduces 
to 


/ da; 

Vdp, 



dW. 


(53) 


One can also arrive directly at (l53l) without defining the 
phase space diffusion (IdlTl by drawing the semiclassical 
limit of Eq. (j49l) and thereby performinga diffusive limit 
analogous to that of a random walker [^. Writing (1551) 
in physical dimensions by using the scales defined in 
Eq. (j4|) finally leads to the Langevin equation of classical 
Brownian motion 1^ 


dp = — 27pdt -I- \/djmkBTenvdW, (54) 


where the dimension of dW is a square root of time. 

It should be remarked here, that the crucial assump¬ 
tion leading to this model is the immediate restoring of 
the pointer state after a jump. The real dynamics will 
take a finite time until the double-peaked wave function 
returns to a pointer state, and there will be several jumps 
in between. We show in the following, by numerically 
simulating trajectories of the pointer state unraveling, 
that this is indeed the case, giving rise to complex dy¬ 
namics. Nonetheless, the general picture of the analytic 
model can be confirmed. 


2. Numerical study 


Using a combination of the Crank-Nicolson method 
and a split-operator technique (ssl - ld^ we generated 
quantum trajectories of the unraveling m and calcu¬ 
lated their first moments x and p and their variances 14 > 
Vp and Cxp- Figure [3] shows exemplarily the temporal 
evolution of p and Vp for two sample trajectories. A sim¬ 
ilar behavior is found in position space (not shown). 

Verification of the stochastic model In order to com¬ 
pare the statistics of the phase space trajectories pre¬ 
dicted by the phase space diffusion (HTl) with that of the 
numerically generated trajectories we consider the expec¬ 
tation values of a finite sample of size N 

1 N 


{Xi 

{Pi 

{xi 


-EJvN)^ (55) 

- VN[p]f , 

- EAr[a:]) (p^ - EAr[p]), 


^n[x] = j^Y1 


N 




Var AT [cc] 
VarAr[p] 

CovAr[a;,p] 


i=l 


I 


N 


N- I 
I 

N-1 

1 

N-1 


E 

2=1 

N 

E 

2=1 

N 

E 
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Figure 3. (a) momentum expectation p and (b) momentum variance Vp of two sample trajectories starting at the same pointer 
state Itt (0, 0)). One observes that in general more than a single jump occurs in the unraveling before the width of the wave 
function is restored to the pointer state width. The calculations are made at k = 50 and all quantities are dimensionless with 
scales defined in Eq. 


where Xi, Pj are the first moments of quantum trajec¬ 
tory i. The sample expectation values (USD are them¬ 
selves stochastic values due to the stochasticity of the 
underlying model. The ensemble variances of Ejv [a;] and 
Eat [p] on the one hand and those of VarAr[a:], VarAr[p], 
and Cov n[x,p] on the other are proportional to 1/iV and 
1/ {N — 1), respectively (see appendix). We compare the 
finite-size variances of the numerically generated trajec¬ 
tories with those predicted by the above model. Figured 
compares Ejv [p\ and Var at [p\ of the sample with the the¬ 
oretically expected deviations. The interval of one stan¬ 
dard deviation around the ensemble mean value Eqo [p] 
should contain Ri 68% of the trajectories of the sample, 
which we confirmed at selected times. The same analysis 
applies for EArjx], VarAr[x], and CovAr[a;,p] and yields the 
same result. Thus, we conclude that the statistics of the 
moments of the pointer state unraveling can indeed be 
described by the stochastic model (1411) . 

Extraction of the diffusion constants Having verified 
the stochastic model we are now able to extract the dif¬ 
fusion constants Dx , Dp and D^p by fitting the variances 
shown in Eqs. (IA5|) to the generated trajectories. Since 
Var[p] only depends on Up, and Cov[a;,p] only on Dp and 
Dxp, fitting is done consecutively by extracting Dp from 
Var[p], Dxp from Cov[a:,p], and finally D^ from Var[x]. 
The numerically obtained values for different k, as well as 
the calculations from the analytic model ((5^ . are shown 
in Fig.[S] We see that the calculated Dp is in good quan¬ 
titative agreement with the simulation. The deviation of 
the calculated D^ and D^p from the simulated ones is 
due to the crude simplification made by the assumption 
that a jump in the unraveling corresponds to a phase- 
space jump of the pointer state without further dynam¬ 
ics; Fig. [3] shows clearly that this is not the case. 

For large n both the analytic model and the numeri¬ 
cal simulation exhibit the classically expected behavior. 
Specifically, the position diffusion D^ and the covariance 
diffusion D^p tend to zero, whereas the momentum dif¬ 


fusion Dp approaches the value 2. 

VI. CONCLUSION 

We identified the pointer states of quantum Brown¬ 
ian motion and derived the stochastic equations of mo¬ 
tion for their trajectories in phase space. The pointer 
states turn out to be Gaussian wave packets with fixed 
widths and a finite position-momentum covariance. The 
phase space trajectory characterizing their motion obeys 
a stochastic differential equation combining momentum 
damping with position and momentum diffusion. In the 
semiclassical limit, it turns into the Langevin equation of 
classical diffusion. This provides us with a consistent and 
transparent picture of the quantum-classical transition as 
induced by quantum Brownian motion. 

We saw that any initial state, which can be represented 
by a superposition of sufficiently separated wave packets 
decays into a mixture of corresponding pointer states, 
and that the ensuing statistical weights are consistent 
with Born’s rule. Using the pointer states and their tra¬ 
jectories one can thus capture an essential part of the 
dynamics of quantum Brownian motion. Our analysis 
depends crucially on linking the nonlinear equation of 
motion for the pointer states to a particular piecewise de¬ 
terministic unraveling of the master equation. It allows 
one to connect the open quantum dynamics to stochastic 
trajectories in phase space. 

The described technique of unraveling a master equa¬ 
tion into pointer state trajectories may become par¬ 
ticularly useful in cases where even approximate solu¬ 
tions of a master equation are not known. For exam¬ 
ple, the full Markovian dynamics of a tracer particle 
in a gaseous environment is described by a complicated 
integro-differential equation not readily accessible by an¬ 
alytic means (2^. If pointer states can be identified for 
this equation we expect their trajectories to be described 
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Figure 4. Eiv [p] and YaxN [p] (solid) as well as their theoretically expected deviations (dashed) due to a finite sample size. The 
fluctuations are well characterized by the dashed lines. Calculations are done at k = 50 for N — 8000 trajectories. 



Figure 5. Dependence of the diffusion constants (solid 
line), Dp (dashed line) and D^p (dotted line) on k. The lines 
are calculated from Eqs. (1S3 of the analytic model, whereas 
the symbols represent the results of the numerical simulation. 
One observes for all three diffusion constants that the ana¬ 
lytic model captures the main features of their dependence 
on K. Specifically, Dp approaches the classical value of 2 in 
the semiclassical regime k —>■ oo, whereas D^ and D^p vanish 
in that limit. 


by a general Chapman-Kolniogorov equation in phase 
space turning into the linear Boltzmann equation in the 
semiclassical limit. 

This work was supported by the DFG via SFB/TR 12. 


Appendix A: Diffusion in phase space 


We analyze the properties of the classical stochastic 
process 


dx = pdt + BiidWi + Bi 2 dW 2 , 
dp = -pdt + B2idWi + B22dW2, 


(Al) 


first, we discuss the ensemble statistics, and afterwards 
the statistics of a sample of finite size. 


1. Ensemble statistics 


In order to derive the statistical properties of the 
stochastic process we recall the stochastic proper¬ 

ties of the real-valued Wiener increments dWi and dIF 2 


E[dW,] = 0, 
dWidWj = Sijdt. 


(A2) 


From these two relations one can derive a differential 
equation for e.g. by expanding to second order in 
dx and keeping only terms to the order df, known as the 
Ito rule of stochastic calculus. One gets 

d (x^) = 2xdx -I- (dx)^ 

= (2xp -b Bf^ + Bf^) dt + 2x (BndWi + Bi 2 dW 2 ) . 

(A3) 

In the same manner one calculates d (p^) and d(xp), 
which then allows the derivation of differential equa¬ 
tions for the ensemble variances Var[x] = E[x^] — E[x]^ 
and Var[p] = E[p^] — E[p]^ and the ensemble covariance 
Cov[x,p] = E[xp] — E[x]E[p]: 

^Var[x] = 2Cov[x,p] -b D^, 

^Var[p] = -2Var[p] -b Dp, (A4) 

^Cov[x,p] = -Cov[x,p] -b Var[p] -b D^^p, 


used in section IV Bl to describe the dynamics of the first 
moments x, p of the pointer states in phase space. At 


where definition (|15)) of the diffusion constants was used. 
Solving the differential equations (IA4I) for initial condi- 
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lions Var[x](0) = Var[p](0) = Cov[a;,p](0) = 0 yields 


Var[a;](t) = + Dp + 2D^p) t - ^Dp (l - e 

— {Dp + 2Dxp) (l ~ e *) , 

Yai[p]{t) = ^Dp{l-e-^*) , 

Coy[x, p\{t) = ^Dp (l - e“*)^ + Dj:p (l - e“‘) . 

(A5) 


can then derive 


d (Ejv[p]^) = 2EAr[p]dEAr[p] + (dEAr[p])^ 


= ( -2E^[p]2 + ^2^ + Bh 


N 


dt 


N 


+ ^Ew[p] ^ (B2idTEi,, + B22dlE2,.), 


2=1 


and finally 


(A7) 


^Var[E^b]] = -2Var[E^[p]] + (A8) 


Here, the ensemble variance of E^v [p] is denoted by 
Var[EAr[p]] = E[Ejv[p]^] — E[EAr[p]]2 and the definition 
of the diffusion constants (H51) was used. 

Similarly, one derives a stochastic differential equation 
for the sample variance of the momentum 


2. Statistics with finite sample size 


We are now interested in the statistics of the first mo¬ 
ments and the variances of a finite sample of size N, and, 
in particular, their deviations around the ensemble val¬ 
ues (jA5ll . With the help of the stochastic model (EH) and 
the sample expectation values (1551) one can write down 
e.g. the stochastic differential equation for the sample 
momentum expectation 


1 . 

dEAr[p] = -EAr[p]dt + ^ X! (-®2idWi,* -f H22dlE2,*), 

^ 2=1 

(A6) 


with independent Wiener increments dWiy and dW 2 y for 
all i. In the same manner as in the previous section one 


dVarAr[p] = {Dp - 2VarAr[p]) dt -f ^ ^ {pi - Eat^) 

X (H2idWi., + H22dW2.,), (A9) 

yielding 

^Var[VarAr[p]] = -4Var[VarAr[p]] -h (AlO) 

for the ensemble variance of VarAr[p]. The variances of 
EAr[a:], VarAr[x] and CovAr[a;,p] are obtained accordingly 
in a long and tedious calculation, giving the solutions 


Var[EAr[a;]] = 

;^Var[a;], 

Var[EAr[p]] = 

;^Var[p], 

Var[VarAr[a;]] = 

2 

N- 

^Var[a;]^ 

Var[VarAr[p]] = 

2 

N- 

-Var[p]2, 

Var[CovAr[a;,p]] = 

1 

N- 

— (Var[a;]Var[p] -|- Cov[a;,p]^) 


(All) 
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